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Ultra-high energy neutrinos are interesting messenger particles since, if detected, they can trans¬ 
mit exclusive information about ultra-high energy processes in the Universe. These particles, with 
energies above 10^® eV, interact very rarely. Therefore, detectors that instrument several gigatons 
of matter are needed to discover them. The ARA detector is currently being constructed at South 
Pole. It is designed to use the Askaryan effect, the emission of radio waves from neutrino-induced 
cascades in the South Pole ice, to detect neutrino interactions at very high energies. With antennas 
distributed among 37 widely-separated stations in the ice, such interactions can be observed in a 
volume of several hundred cubic kilometers. Currently 3 deep ARA stations are deployed in the ice, 
of which two have been taking data since the beginning of 2013. In this publication, the ARA detec¬ 
tor “as-built” and calibrations are described. Data reduction methods used to distinguish the rare 
radio signals from overwhelming backgrounds of thermal and anthropogenic origin are presented. 

Using data from only two stations over a short exposure time of 10 months, a neutrino flux limit 
of 1.5 X 10“® GeV/cm®/s/sr is calculated for a particle energy of 10^® eV, which offers promise for 
the full ARA detector. 


I. INTRODUCTION 

In 1966 Greisen, Zatsepin and Kuzmin predicted an 
interaction of ultra-high energy cosmic rays (UHEGRs) 
with the recently discovered cosmic microwave back¬ 
ground radiation mm- In such interactions, pions are 
produced resonantly, subsequently decaying into neutri- 
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nos, as first postulated by Berezinsky and Zatsepin in 
1968 [3]. Due to the contribution of the Delta resonance 
to the cross section for this interaction, UHECRs are un¬ 
able to reach us from sources on cosmological distance 
scales, i. e. beyond tens of Mpc. This implies a sharp 
cutoff in the cosmic ray spectrum at an energy of around 
IQiS-® eV, which has been confirmed by the largest cos¬ 
mic ray air shower detectors, Telescope Array H] and 
the Pierre Auger Observatory [5]. However, the lack of 
UHECRs arriving at Earth could be due to a number of 
underlying reasons dependent on their properties. The 
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mass composition, for example, plays a crucial role in 
determining the dominant energy loss for UHECRs [6]. 
Furthermore, this energy loss is influenced by the distri¬ 
bution of sources and the primary energy spectrum of 
cosmic rays. 

As a consequence, the neutrino flux depends strongly 
on all three parameters: the UHECR composition, their 
energy spectrum and their source distribution. Hence, a 
measurement of this flux can be used to place constraints 
on those parameters. Moreover, due to the Greisen- 
Zatsepin-Kuzmin (GZK) effect and similar absorption 
mechanisms for HE gamma rays, neutrinos are the only 
feasible known particles for the study of UHE sources 
more distant than a few tens of Mpc. Neutrinos, charge¬ 
less and only weakly interacting with extremely low cross 
sections, arrive at detectors unscattered and undeflected 
by intervening particles and fields, and may be correlated 
to UHE sources at the furthest distances in the Universe. 

Many other UHE neutrino source models have been 
offered. Recent IceCube results challenge some of these 
models that predict larger neutrino fluxes at EeV ener¬ 
gies, given the somewhat smaller-than-predicted number 
of neutrinos measured by IceCube with energies in the 
PeV range [7]. Nevertheless, other theoretical models 
predict EeV neutrino fluxes that comply with IceCube 
limits, e.g. from pulsars [5], blazars [5] or from the after¬ 
glow of gamma-ray bursts [TOl [H]. A more detailed dis¬ 
cussion of such theoretical models will not be presented 
in the framework of this experimental paper, but may be 
found in the references. 

The expected flux of GZK neutrinos at Earth from dif¬ 
ferent cosmic ray models is very low m and, in combi¬ 
nation with the low interaction cross section m, leads to 
an interaction rate of less than 1 GZK neutrino per giga- 
ton of matter per year. Therefore, large detectors, cover¬ 
ing several hundred cubic kilometers of water equivalent 
matter are needed to record neutrino events in sufficient 
quantity to investigate their flux. 

The large attenuation length of Antarctic ice to ra¬ 
diofrequency waves, of 0{1 km), opens the possibility to 
space detectors on a comparable scale and to utilize co¬ 
herent radio emission from neutrino induced cascades in 
radio transparent media; the so-called Askaryan effect 
[nms], which has been verified in various experiments 
[MU]- In the interactions of high energy neutrinos with 
electrons or nuclear matter, electromagnetic (EM) cas¬ 
cades are produced which build up a net negative charge 
of roughly 20% close to the shower maximum. This im¬ 
balance originates mainly from Compton scattering of 
cascade photons on atomic electrons. Smaller contribu¬ 
tions are added by other ionizing effects such as positron 
annihilation with atomic electrons [HI [20] ■ The net 
charge acts as a moving current and emits electromag¬ 
netic waves, which become coherent at wavelengths com¬ 
parable to the lateral cascade dimensions. This is valid 
in the radio regime. In the case of coherent emission, the 
strength of the EM far field is proportional to the cascade 
energy. The frequency spectrum of the Askaryan signal 


depends strongly on the observation angle. As described 
in [19], the strongest signal is observed at the Gherenkov 
angle for frequencies around 1 GHz. The signal distribu¬ 
tion around this angle can be approximated by a narrow 
Gaussian distribution. At lower frequencies the signal 
is weaker but the angular distribution is much broader. 
This trade-off needs to be considered when designing a 
neutrino detector utilizing this emission. 

Askaryan radiation is a consequence of a neutrino of 
any flavor interacting in a radio transparent medium such 
as ice. However, as the effect ultimately comes from the 
induced EM cascade, the detection efficiency is strongly 
dependent on neutrino flavor. Gharged-current inter¬ 
actions and Ve+e~ elastic scatters convert a large fraction 
of the neutrino energy into the EM cascade. Neutral- 
current interactions of all v flavors may initiate hadronic 
cascades which receive, on average, only 20% of the neu¬ 
trino energy. These hadronic cascades themselves will 
produce EM sub-cascades due principally to decaying 7r° 
particles, but a fraction of energy is lost to hadrons. On 
the other hand very high energy EM cascades undergo 
elongation via the Landau-Pomeranchuk-Midgal (LPM) 
effect Emm which does not affect hadronic cascades. 
Radiation emitted by the LPM-elongated EM cascades 
is strongly beamed at the Cherenkov angle and thus is 
less likely to intersect detection antennas. Therefore, de¬ 
spite their stronger signal emission, Ve are not necessarily 
the dominant detected flavor given an astrophysical ex¬ 
pectation of (zze : : Vt) = (1:1:1) and most of the 

acceptance is expected to arise from hadronic cascades. 


II. THE INSTRUMENT 

The Askaryan Radio Array (ARA) is a neutrino detec¬ 
tor, currently under construction at the geographic South 
Pole next to the IceCube experiment. It is designed to 
utilize the Askaryan effect to detect interactions of GZK 
neutrinos in the South Polar ice sheet. At the chosen site, 
thousands of square kilometers of ice with a thickness of 
about 3 km are available to act as a radio transparent 
detector medium and to allow for the construction of a 
0(100 Gt) detector. Due to its low temperature, between 
—55 °C and — 30 °C in the top 2 km [23], the South Pole 
ice sheet has low radio attenuation. On average, an at¬ 
tenuation length of 820 m integrated over the top 2 km 
of ice has been measured for frequencies around 300 MHz 
[24] . Furthermore, the Amundsen-Scott station provides 
the infrastructure to support large projects such as the 
ARA experiment. 


A. General design 

The ARA detector baseline consists of 37 antenna clus¬ 
ters ^‘stations”) spaced by 2 km in a hexagonal grid (Fig¬ 
ure [^. Each station is designed to operate as an au¬ 
tonomous neutrino detector and simulations have shown 
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Figure 1. An area map of the planned ARA detector at South 
Pole. The stations are indicated by the black circles. Red 
filled circles denote the currently deployed stations. 
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a grid spacing of 2 km to nearly maximize the array’s 
effective area at neutrino energies of 10^® eV |24j . 

Each station comprises 16 measurement antennas, de¬ 
ployed on strings in groups of 4 at the bottom of 200 m 
deep holes. In the baseline design, these antennas form 
a 20 m X 20 m X 20 m cube (Figure |^. This design is 
in the process of being optimized based on analysis re¬ 
sults from the first ARA stations and simulations. Each 
hole contains two antennas of horizontal and two anten¬ 
nas of vertical polarization, all recording data between 
150 MHz and 850 MHz. Two separated polarizations are 
chosen to be able to determine the polarization of the 
incoming signal, which is important for neutrino recon¬ 
struction. The antenna names are composed of the string 
number as -D#, their position on the string (T for top, 
B for bottom) and their polarization (V for vertical, ff 
for horizontal). 

The antennas are deployed at depths between 170 m 
and 190 m to minimize the effects of ray-tracing in the 
ice. Due to the depth-dependence of the temperature and 
density of the South Pole ice sheet, the index of refraction 
changes with depth [25]. This effect is strongest in the 
top 200 m, starting from an index of 1.35 at the surface 
and changing to a value of 1.78 for the deep ice at a depth 
of roughly 200 m. As described in [211 US] , this causes 
the path of radio rays to be bent downwards rendering 
vertex reconstructions difficult. Moreover, a shadowed 
area is produced where signals cannot reach shallowly 
deployed antennas, thus reducing effective neutrino vol¬ 
ume. Therefore, a deep deployment of the antennas is 
favorable. 

In addition to the receiver channels, 4 calibration 
transmitter antennas are deployed on two extra strings, 
D5 and D6. These pulsers can transmit transient signals 
or continuous broadband noise for calibration of the sta¬ 
tion timing, geometry, and signal efficiency. Transient 
emissions are tied to a GPS clock which allows separa¬ 
tion of them from other recorded RE signals by timing. 
The pulsers are positioned at a distance of approximately 
40 m to the station core at a similar depth as the mea- 


Figure 2. The baseline design of an ARA station with a zoom 
illustrating the string details and a view of the deployed an¬ 
tennas of both polarization. 


surement antennas. Each hole contains one antenna of 
each polarization. 

The antennas used in ARA are birdcage dipoles for the 
vertical polarization (Vpol) and ferrite loaded quad-slot 
antennas for the horizontal polarization (Hpol). Given 
the drilled antenna holes with a diameter of only 15 cm, 
the design of Hpol antennas with reasonable sensitiv¬ 
ity down to 150 MHz is very challenging. Slotted cop¬ 
per cylinders show reasonable low-frequency performance 
with a voltage standing wave ratio below 3, for frequen¬ 
cies above 300 MHz [ 22 ] • This can be further improved 
by adding ferrite material in the cylinder core. 

The signals recorded by the antennas are first filtered 
by a bandpass and notch filter, to reject frequencies out 
of band as well as narrow-band communications. After 
filtering, signals are amplified by Low Noise Amplifiers 
(LNAs) and transmitted analog to the surfa ce th rough 
fiber cables via optical Zonu links (Figure 3(a)). At 
the surface, signals are filtered again, split and fed to 
the trig ger sy stem as well as to the digitization system 
(Figure r3(b)[ ). A first calibration of the full signal chain 
and antenna response after deployment in the ice can be 
found in Section III PI 

In the trigger system, the signal is processed by an 
integrating tunnel diode, producing energy envelopes of 
the incoming waveforms, which can be processed in the 
trigger electronics mounted on the Triggering Daughter 
board for ARA (TDA). On this board, the signal is read 
into the digital electronics and logic, implemented in an 
FPGA. This logic then determines whether the event sat¬ 
isfies the trigger condition. This condition is currently a 
simple multiplicity trigger, requiring signal on 3 out of 
8 channels of one polarization. Investigations are being 
performed to replace this with a smarter algorithm which 
provides more efficient background rejection and better 
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Figure 3. Left: The components of the down-hole signal chain 
on each string in the ARA stations. Right: The surface Data 
acquisition system of the ARA stations, showing the most 
important components. Components framed in yellow are 
common to all strings. 


Figure 4. The cumulative live time of stations A2 and A3 in 
2013. Horizontal line segments indicate extended downtimes 
of the detectors. 


B. Performance of the ARA stations A2 and A3 


signal retention. 

The digitization system is located on the Digitizing 
Daughter board for ARA (DDA). In this system the data 
are sampled by the IRS2 ASIC, a digitization chip capa¬ 
ble of sampling data at a rate of 4 GS/s. In the ARA 
detector, the sampling speed is tuned to 3.2 GS/s. The 
IRS2 chip contains 8 channels each with a 32k-element 
Switched Capacitor Array (SCA). The 32k elements are 
further subdivided into 512 randomly write-addressable 
blocks of 64 samples each. Analog sampling is continuous 
and is stopped by an external trigger to signal the start 
of digitization and readout of the analog storage blocks 
of interest. Performance of this early version of SCAs 
with deep analog storage buffers shows promise for mul¬ 
tichannel high-speed, low-power samplers. Power con¬ 
sumption is in the range of 20 mW per channel. In prin¬ 
ciple deadtimeless operation is possible due to the deep 
analog buffer, however, noise issues related to simulta¬ 
neous readout and digitization have prevented operation 
in this mode to date. The calibration of this digitizer is 
presented in Appendix [A} 

Both the TDA and DDA boards are mounted on the 
ARA Triggering and Readout Interface (ATRI) which 
provides all logic for the data acquisition systems in a 
single Spartan-6 FPGA. This FPGA is programmed by, 
and exchanges its data with, an Intel atom-based Single 
Board Gomputer (SBG) which handles the data trans¬ 
fer to storage on disks in the IceCube Laboratory (IGL). 
Currently an event rate of 5 Hz can be accommodated on 
a USB link between the FPGA and the SBG. In Decem¬ 
ber 2015 a much higher bandwidth PGI Express bus was 
installed to increase the acceptable station trigger rate. 

The ARA detector in its current form consists of 3 
stations, of which two were taking data in the year 2013 
and produced the data for this analysis. 


The first data from the deep ARA detector have been 
recorded by stations A2 and A3. Station A1 could not 
deliver data in the year 2013 due to an issue in the com¬ 
munications system. This problem was repaired in the 
2015-2016 Antarctic summer season and all 3 stations 
are now operational. 

The positions of A2 and A3, as embedded in the full 
ARA37 design, are shown in Figure Their structure 
closely follows the ARA baseline design as described in 
Section p^I A[ After being deployed in February 2013, the 
stations recorded data for 10 months until the end of that 
year (Figure]^. Due to various infrastructural issues and 
optimizations which interrupted the detector operation, 
there were several extended periods of down time, some¬ 
times lasting for days. Therefore, the two detectors were 
only running 75% of the time and correspondingly accu¬ 
mulated about 228 (A2) and 220 (A3) days of live time 
during those months. Meanwhile, station operation has 
become more stable through debugging and optimization 
of the data acquisition (DAQ) firmware and software. In 
addition to that, new monitoring tools have been devel¬ 
oped which allow us to identify and solve problems within 
a few hours. This resulted in a significant rise in live 
time for the year 2014. The dead time during operations 
due to digitizer occupancy and the limited data transfer 
bandwidth is very small and less than 1% of the total 
run-time. 

Of the 32 deep in-ice measurement channels in the two 
stations, 31 are fully operational. The bottom Hpol chan¬ 
nel on string 4 (D4BH) in A2 shows strong noise fluctu¬ 
ations which are believed to be due to a damaged LNA. 
Figure shows pulser waveforms of both polarizations as 
recorded by A3. One can see the clear absence of a sig¬ 
nal in the antennas polarized perpendicular to the emit¬ 
ted signal, which indicates that the two polarizations are 
very well-separated. The RMS of the background noise 
on the two stations is relatively stable throughout the 
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Figure 5. Signals emitted by a Vpol (a) and an Hpol (b) 
pulser, as recorded by station A3 string D3. As is evident, 
the polarization separation is very clean and only the antennas 
of the emitted polarization show a response. 


year, as shown for A2 in Figure 



Figure 6. The RMS of events recorded by 4 selected mea¬ 
surement channels from A2 in 2013, plotted in mV versus the 
day of 2013. The large variation observed in channel D4BH 
is indicative of failure of that particular channel. 
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C. Calibration of station geometry and timing 

After having achieved a stable timing from the digi¬ 
tizer chip, the systematic errors in this timing and the 
precise positions of the antennas in a station need to 
be determined, to allow for accurate vertex reconstruc¬ 
tions. The antennas are suspended on four strings in 
four vertical holes, connected by stiff cables. Their XY- 
coordinates can thus be assumed to coincide within one 
hole. Furthermore, the vertical distances and the cable 
delays have been measured and are assumed to be cor¬ 
rect with a negligible error. Parameters which still need 
to be calibrated are the position of each string and the 
relative time delay between them. One string has to be 
chosen as perfectly positioned and the rotation of the 
station around this string needs to be fixed to obtain a 
well determined coordinate system. With these assump¬ 
tions there remain 17 parameters of positions and cable 
delays to be calibrated (neglecting uncertainties in the 
index-of-refraction model). 

Such calibration is performed by using calibration 


Figure 7. Results of the fits for geometrical and timing cali¬ 
bration of the strings in station A3. Note that string D3 has 
been used as the reference and is fixed in the -minimization 
of this calibration. 

pulser signals and determining the arrival time difference 
between signals on different antennas. For this quantity 
and the geometrical positions of the antennas one can set 
up an equation for each possible antenna pair and for all 
4 pulsers. Considering the two polarizations, a system of 
28 independent equations can be constructed. For these 
equations we can set up a value as 

='^[ - dt%,rei) 

-\-X]^ * “h * ‘^yij Z}^ * 2^2^ 

-tfe,ref ■ 2c^dtk^J - rf + rl ] . (1) 

Here we use c = 0.3/1.755 m/ns as the speed of light 
in ice at the average antenna depth, taken from [25] . The 
coordinates of the pulsers are denoted by x, y, z and the 
arrival time difference by dt, with the subscripts k for 
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Figure 8. The residual of the rooftop pulser reconstruction 
with A3 (a) before and (b) after the geometrical calibration. 
The bimodal distributions contain noise events with a resid¬ 
ual of roughly ® and signal with lower residuals which 
migrate to smaller values after calibration. 


the used pulser and j for the respective antennas. The 
parameter r indicates the distance of the antenna to the 
station center. The presented is closely related to 
the equations which are used for the reconstruction al¬ 
gorithm, described in Section IV B| Standard minimizer 
tools are used to minimize the using multiple and 
different seed values for the input parameters. The av¬ 
erage of the outcomes is taken as the final result. Fig¬ 
ure shows how different parameters are constrained by 
this method for station A3. It becomes apparent that 
the current geometrical setup is strongest in determining 
the string depth while the constraints on the X- and Y- 
position, and the inter-string delays are relatively weak. 


The result of this calibration is checked via the recon¬ 
struction of an independent pulsing antenna, mounted at 
a distance of about 4 km from both ARA stations on the 
rooftop of the IceCube Laboratory. The reconstruction 


algorithm described in Section IV B is used for this cross 


Figure 9. The directional rooftop pulser reconstruction with 
A3 (a) before and (b) after the geometrical calibration. The 
axes show the difference between the reconstructed and the 
true azimuth (X-axis) and zenith (Y-axis) angle. Reconstruc¬ 
tion quality criteria are only applied loosely. It should be 
noted that the data from the rooftop pulser are not part of 
the calibration data sample. 


check. The figures of merit are the stability of the re¬ 
construction as well as the residual which indicates the 
internal consistency of the station geometry. Signals from 
the rooftop pulser are not tied to a GPS clock and have 
to be filtered out of all recorded data by other means. 
The residual is plotted before and after calibration for 
all data recorded during a rooftop pulser run in Figure 
Two peaks are visible in this distribution: one for 
noise waveforms with high residual around 10“^'® and 
one for signal which shifts to significantly lower values 
after the calibration. This indicates that the assumed 
geometry is more consistent with the measured timing 
after the calibration has been applied. Figure shows 
the actual result of the reconstruction compared to the 
expected value. After calibration, the reconstruction is 
much more self-consistent than before. Currently, in the 
two stations all operating channels but two in A2 are 
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calibrated in position and timing in this way. The two 
omitted channels in station A2, D3BV and D3BH, show 
a puzzling timing offset of several nanoseconds which is 
not corrected in the presented calibration. The source of 
this offset is unknown to date, but more measurements 
have been taken to further improve the precision of this 
calibration. 


D. Signal chain calibration 

In the amplitude calibration of the ARA signal chain, 
we try to determine the noise temperature of the en¬ 
vironment, the noise figure of the signal chain and the 
directional gain of the antenna. 


1. Determination of ambient noise 

The ARA antennas are exposed to various sources of 
radio noise. The main contributors are the ice surround¬ 
ing the antenna, with a depth dependent temperature 
between 220 K and 270 K m and sky sources like the 
atmosphere and the Galactic center. The contributions 
of the Sun, moon and cosmic background radiation are 
negligible for the ARA antennas. Also the contribution 
from the bedrock under the ice, with a temperature of 
roughly 273 K, could be determined to be negligible due 
to the radio attenuation in ice and the disfavored incom¬ 
ing angle in the directional gain pattern of the antennas. 

The radiation from the atmosphere and galactic noise 
both approach the ice from the top and enter the ice ac¬ 
cording to Snell’s law at the boundary of a medium. The 
minimum elevation angle, as viewed from the in-ice an¬ 
tennas, above which sources emanating from above the 
ice can be observed by those in-ice antennas is approxi¬ 
mately 55°. Atmospheric and galactic temperature pro¬ 
files have been extracted from PH HH] ■ The attenuation 
in the atmosphere is assumed to be negligible and the at¬ 
tenuation through the top ice is normally low. However, 
the contribution to the noise of the antennas is relatively 
small, with an equivalent temperature of Tgky = 18.3 K 
at 300 MHz for Vpol antennas. This is mainly due to the 
steep incoming angle which is highly disfavored by the 
directional gain pattern of the antennas. 

To calculate the power spectrum received from the ice 
one has to consider it as divided into semi-transparent 
volume-elements. The brightness B' of each element is 
the brightness of a black body B, reduced by the limited 
emissivity e of the ice as: 

B'{u)=e-B{u) = --B{v), (2) 

a 

with a being the attenuation length at 300 MHz at the 
given depth. This equality is valid under the assumption 
of thermodynamic equilibrium in the ice |29j . The total 
noise power at the antenna is the integral over all volume 


elements in a given range i? around the antenna 

f R 

^27rA^ / j B'{v)ij{r,d)G{y, 9) cos 9dOdr 
Pioe{y) = --, (3) 

27r f G(iy, 9) ■ cos 9d9 

~ 2 ' 

with the wavelength A, the directional antenna gain G 
and the reduction factor rj due to the finite attenuation 
length for radio waves in ice. The factor 1/2 is applied to 
account for the single polarization of the antenna. In the 
integral of Equation the contribution of radiation re¬ 
flected at the surface has to be taken into account, which 
accounts for roughly Tren = 53 K of noise power. 

The total noise temperature seen by an ARA antenna 
can be calculated to be 

Tant = Pice + PreB + Psky (4) 

This results in 247 ± 13 K for Vpol and 249 ± 13 K for 
Hpol antennas at 300 MHz. The main contributions to 
the errors come from uncertainties in the ice tempera¬ 
ture measurement, the critical angle for surface reflection 
and the antenna gain pattern. The temperature value 
currently used in the ARA simulation is calculated for 
one single frequency of 300 MHz. Since this temperature 
changes over the ARA frequency range by up to 30 K, 
the latter will be assumed as a systematic error for fur¬ 
ther calculations rather than the smaller calculated error 
of 13 K. 


S. Signal ehain ealibration 

The calibration of the signal chain is particularly chal¬ 
lenging since the deployment of fixed strings in 200 m 
deep holes allows neither a detailed investigation of the 
directional gain pattern of the antenna nor noise figure 
measurements in-situ. 

The noise factor F of the signal chain has been mea¬ 
sured for each channel at the surface at room temper¬ 
ature prior to deployment. This noise factor is used to 
calculate the noise temperature of the signal chain. The 
change due to the signal chain being lowered into the ice 
with a local environmental temperature of roughly 220 K 
at a depth of 180 m is estimated as a linear change with 
ambient temperature. The total noise contribution de¬ 
pends on the transmission coefficient t of the antenna. 
For perfect coupling it can be calculated to be (e.g., at 
300 MHz): 

Ttot =Ara„t+r(180 m).(E-l) 

= 247K-h220 K-(1.6 - 1) = 379 K, (5) 

where t = 1 for perfect coupling. 

With knowledge of this noise floor, the directional gain 
can be determined using the external noise sources at 
each station. These sources are noise diodes, emitting 
a flat power spectrum in the relevant frequency range, 





Figure 10. Power plotted versus frequency for a recorded 
Vpol signal from an external source (blue solid line), with a 
fiat input spectrum between 150 and 1000 MHz. The cumu¬ 
lative power distribution shows that most of the signal power 
is recorded below 500 MHz. The green dashed line shows the 
recorded spectrum for the natural ambient noise without ad¬ 
ditional signal for comparison. 


which can be attenuated by up to 30 dB. They are con¬ 
nected to the pulser antennas, installed in the vicinity of 
the ARA stations. Figure shows the recorded power 
spectrum for a typical receiver antenna and a pulser with 
a non-attenuated input source. This figure further illus¬ 
trates the sensitivity of the signal chain in different parts 
of the frequency spectrum via a cumulative power distri¬ 
bution. Most signal power is recorded below 500 MHz in 
this measurement. 


tor from 

47rR / Pout - 

X V UUPtNout 

■ \/ir-A^ant + 220 K-(F-l), (6) 

with A"ant being the power corresponding to Tant, tr/t 
the transmission coefficient between a given antenna and 
the signal chain (r for receiver, t for transmitter), Pt the 
input power to the noise source antenna and R the dis¬ 
tance between transmitter and receiver antennas. A^out 
is the measured output power without any applied signal 
and Pout denotes the total recorded power when the noise 
sources are operating. These are the experimentally mea¬ 
sured values. The transmission coefficient t is taken from 
XFDTD simulations [3T] of the deployed Vpol antennas 
in ice and from NEC2 simulations of the Hpol antennas. 
For the directivity, we assume for now = Dt = D, for 
a transmitter and receiver of the same polarization. 

The antenna gain can then be calculated as 



G = D-t. 


( 7 ) 


Figure [T^ shows the gain measured for the A3 Vpol 
antennas and the A2 Hpol antennas, both plotted with 
systematic errors, which are dominant over statistical 
uncertainties in this measurement. They are derived 
from uncertainties in the input spectrum, the ambient 
noise power and the used noise figures. The comparison 
to simulations shows that differences, especially beyond 
500 MHz, are not covered by our current understanding 
of these uncertainties. This may be due to an imperfect 
simulation of the antennas in ice, which is indeed very 
challenging, or due to unaccounted error sources in the 
measurement. The signal strength at frequencies beyond 
500 MHz is relatively low (see Figure 10), which is ex¬ 


pected due to the lower sensitivity of the signal chain 
at high frequencies. Such behavior allows for a stronger 
influence of non-linear effects in this region which are 


Under the assumption that all installed antennas of 
the same polarization have the same angular gain pat¬ 
tern, the geometric relation between the pulsers and all 
8 measurement antennas of one polarization in both sta¬ 
tions can be used to measure the gain pattern for different 
angles. For this measurement, we assume further that the 
directionality in azimuth is isotropic. With the available 
antennas we can establish 24 different relations in eleva¬ 
tion which are clustered at distinct angles, as shown by 
the example of A3 in Figure El In fact, for station A3 
all bottom antennas are at the same depth as the two 
pulsers (D5 not shown), while the top antennas are ele¬ 
vated by an angle of roughly 25°. In station A2, the one 
fully operational pulser is mounted at a similar depth to 
the top antennas, while the bottom antennas receive the 
pulser signal at an inclination of —25°. The antenna di¬ 
rectivity D at a given angle is then calculated via the 
Friis transmission equation m and the known noise fac- 



Figure 11. The relative geometry of the Vpol antennas to the 
D6 pulser in A3. 
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derived for a Vpol antenna in ice based on an adapted 
NEC2 simulation. 

The visible asymmetry in zenith angle, strongest at 
350 and 400 MHz, is expected to be due to the connect¬ 
ing cable which is fed vertically through the antennas to 
a connector in the center. This effect has been observed 
in calibration measurements in air. Its strength in ice 
could not be quantified in the presented measurement 
since antennas of one type (bottom/top) have not been 
measured at positive and negative angles. Furthermore, 
a difference in the source could influence the current pic¬ 
ture since certain angular combinations only appear for 
a given source. 

The wide clustering of data points is most likely caused 
by the change of input source and the difference between 
top and bottom antennas. Furthermore, the shape of the 
hole walls around the antennas may have a signihcant 
influence which still needs to be quantified. 


3. Discussion 


Figure 12. Sample measurements of the ant enna’s directional 
gain at the available angles (see Figure 111 for the two ARA 
stations and polarizations: (a) A3 Vpol, (b) A2 Hpol. The 
data are plotted on a linear scale versus frequency. For com¬ 
parison the current status of the simulation is shown as a 
black solid line. String D4 on A3 was not operating during 
the time of the measurement. The black dashed line repre¬ 
sents the lower limit on the signal gain used to derive the 
systematic error on the detector sensitivity. 


difficult to quantify. One possible source for such non- 
linearities is the digitizer chip and its calibration. In Ap¬ 
pendix it is shown that non-linearities are observed 
in the ADC-to-voltage conversion gain. In addition, the 
sample timing can have errors of 0(100 ps). Such imper¬ 
fections may add artificial power to the Fourier spectrum 
depending on the level of disturbance, as illustrated in 
Figure This could explain the excess in the antenna 
gain measurement at frequencies above 500 MHz. The 
influence of non-linearities on the frequency range below 
500 MHz, given sufficient signal strength, is however ex¬ 
pected to be very small. The visible difference between 
top and bottom antennas is due to cable feed-through. In 
the ARA stations, an antenna contains a vertical cable 
feed-through for each antenna which is mounted below. 
This has a significant influence on the angular gain pat¬ 
tern, which has not yet been fully simulated. 


Figure [T^ shows the gain measured by all possible Vpol 
antennas versus the reception angle. As is evident, the 
antennas are not distributed equally over the angular 
range, but concentrated at a few points. Therefore, it 
is difficult to make a prediction for the full angular re¬ 
sponse. The gain pattern currently used for the bottom 
antennas in the ARA simulation (Section |HI[ ) is included 
in Figure 13 for comparison. This simulation has been 


The calibration presented above is an initial step to un¬ 
derstand the behavior of the ARA signal chain and the 
ambient noise. Several discrepancies between the mea¬ 
surement and the currently used simulation have been 
pointed out which exceed the shown systematic errors 
on the measurement. These discrepancies need to be re¬ 
solved, to obtain a comprehensible frequency spectrum of 
recorded signals. Plans to improve the shown measure¬ 
ments are currently under development. One possibility 
to achieve a better understanding of the antenna gain is 



Figure 13. The directional gain results versus reception an¬ 
gle from A2 and A3 for all Vpol antennas compared to the 
current simulation of the bottom antennas (green line) at dif¬ 
ferent frequencies. All data are normalized to an isotropic 
directionality pattern. The three data points per antenna 
originate from the three available calibration sources used in 
this calibration. Gains around —25° angles can only be mea¬ 
sured in the A2 geometry, gains around -|-25° only in A3. 
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to move pulsers vertically in the hole and to take mea¬ 
surements at different depths. 

From the signal power distribution in this measure¬ 
ment we conclude that the calibration values are only 
reliable in a frequency range of up to roughly 500 MHz. 
Above that frequency we consider our signal chain to be 
understood more poorly. 

In the lower part of the spectrum, up to roughly 
500 MHz, simulations appear to underestimate the Vpol 
antenna gain at the measured angles and to overestimate 
Hpol antennas. An average underestimation for Vpols of 
15% and an overestimation for Hpol antennas of 30% can 
be obtained between simulation and measurement. Given 
this, we choose the negative systematic error on the an¬ 
tenna gain to be 0 for the Vpol antennas and —30% for 
the Hpol antennas. 

To account for all errors determined in this signal chain 
calibration we combine them into an uncertainty on the 
signal to noise ratio (SNR), which is used to determine 
the error on the detector sensitivity (see Section VII. The 
signal to noise ratio in power can be calculated to be 


SNR = 


t-D-R 


Slg 


AN, 




( 8 ) 


with the transmission coefficient t, the directivity D, the 
incoming neutrino signal Psig, the ambient noise 
and the signal chain noise Ngc- The uncertainty on the 
signal chain noise figure has been measured to be 10% on 
average, and the general uncertainty on the transmission 
coefficient is assumed to be 10%. As mentioned in section 
the error on the ambient noise temperature is taken 
30 K. All error values are summarized in Table HI 
Given those values the resulting relative error on the SNR 
in power is —32% for Hpols and —10% for Vpols. 


IID1 
to be 


Table I. The estimated errors from various sources included 
in the lower systematic error on the SNR in power. 

Source estimated error 

Directivity D 30% (Hpol), 0% (Vpol) 

Transmission coefficient t 10% 

Signal Chain N^c ~ 10% (frequency dependent) 
Ambient noise temp Aant 30 K 

Total_ S2%{Hpol),10%{Vpol) 


Since we do not have conclusive calibration results for 
the upper part of the spectrum, we assume a worst case 
scenario for the lower systematic error and apply a low 
pass filter on the signal in this area. To quote an upper 
systematic error, the error sources in the upper frequency 
range need greater understanding than has been achieved 
at this point. Therefore, no upper systematic error is 
presented for now. This, however, has no impact on the 
flux limits presented later in this document. 

The relative lower limit on the SNR is represented by 
the dashed black line in Figure 12 For the calculation 


of the systematic error on the detector sensitivity, this 
lower limit gain is applied to simulated neutrino signals 
to estimate its influence on the effective area (see Section 

Vl. 


III. SIMULATIONS 

The simulation of neutrino vertices for ARA is per¬ 
formed with the AraSim code, which is described in de¬ 
tail in |26l [32]. In this section a short summary of the 
simulation is presented. 

In AraSim, forced neutrino interactions are generated 
uniformly over a cylindrical volume. This volume is cen¬ 
tered on the simulated ARA station with incoming neu¬ 
trino directions uniform in cos 6*. It is bounded by the 
bedrock under the ice in depth and by an energy de¬ 
pendent radius, chosen to include all possibly triggering 
events. For each interaction, a weight is calculated based 
on the probability that the neutrino would interact at 
the given point after having passed earthbound material 
along its trajectory. This probability depends on the en¬ 
ergy dependent cross-section for a neutrino interaction 
[33] and the summed number of radiation lengths along 
the neutrino path. The primary neutrino energy spec¬ 
trum can be chosen freely. 

For each simulated neutrino interaction, a cascade and 
its radio frequency emission are modeled from theoretical 
approximations. The results obtained in the following 
analysis are based on modeling of the Askaryan emis¬ 
sion in the frequency domain according to [20j . which 
has since been updated to a semi-analytical method of 
simulating the emission for each event based on |34] . 

The trajectory of the RF signal to the antenna is cal¬ 
culated using fitted models for the index of refraction 
pS] and the final signal strength is derived taking into 
account the depth dependent attenuation length 



Figure 14. The effective area of the two ARA stations as a 
function of neutrino energy. 
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Finally, the signal chain and the trigger system are 
modeled after calibration measurements of their compo¬ 
nents in the laboratory. For the antenna response, a 
NEC2 simulation is used. The trigger logic is modeled 
in AraSim as it is currently set up in the ARA stations: 
whenever 3 out of 8 antennas of the same polarization 
cross a given threshold, an event is recorded. 

Thermal background noise is modeled using the av¬ 
erage frequency spectrum from unbiased forced trigger 
events recorded throughout 2013, measured with each an¬ 
tenna in station A3. The production of a significant sam¬ 
ple of thermal noise data is difficult due to the 6 ct power- 
threshold which is currently used in the trigger system to 
limit the event rate to 5 Hz. Therefore, simulated noise is 
exclusively used to develop and initially test algorithms 
while final checks and the estimation of background are 
performed on recorded data. 

From simulations, the effective area of the two ARA 
stations can be calculated at the trigger level to be: 


A,s{E) 


VgenjE) 1 


(9) 


where V^en is the cylindrical volume over which events are 
generated, A^gen is the number of events which have been 
generated, uji is the weight of each event i which triggered 
the detector and Lint is the interaction length at the given 
energy. The effective area for the combined A2 and A3 
detector is plotted at the trigger level in Figure In 
this plot, the single effective areas per station are also 
shown as well as the effective area for events which are 
coincident to both stations. These coincidences amount 
to roughly 5% of all events at an energy of 10^® eV. 

The simulated dataset used in the present analysis con¬ 
tains equal numbers of neutrinos in quarter-decade en¬ 
ergy bins between 10^® eV and 10^^ eV. All figures, show¬ 
ing simulated neutrino events are based on this sample. 


IV. DATA ANALYSIS 
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Figure 15. (a) The hit pattern for a simulated event contain¬ 
ing a signal waveform. The hits are indicated by the colored 
squares. Hits which correspond to the signal wavefront are 
marked by black dashed ellipses to underline their visible reg¬ 
ular pattern, (b) A typical hit pattern for pure thermal noise. 


The analysis of the first data from A2 and A3, recorded 
in 2013, has been optimized for sensitivity to neutrino 
interactions at a fixed rejection of thermal and anthro¬ 
pogenic backgrounds. 

The ARA detector records events at a rate of roughly 
5 Hz. These events are mostly thermal noise and to a 
lesser degree, backgrounds of anthropogenic origin. In 
this analysis these backgrounds are reduced in two steps: 
with a thermal noise filter and by application of angu¬ 
lar cuts to reconstructed vertices. All algorithms have 
been developed and tested on a 10% subset of the full 
recorded data, in the following called “burn-sample”, to 
avoid a bias in the analysis. Cuts are developed to re¬ 
duce the expected background to approximately a factor 
10 beneath the level of expected neutrino events. After 
the cuts are finalized, the analysis is applied to the full 
recorded data set for the year 2013. 


A. Thermal noise filtering 

Thermal noise filtering is performed with the time se¬ 
quence filter, developed for the close-to-cubical ARA sta¬ 
tion geometry. Further details about a first version of the 
method, described below, can be found in |35) . 

The algorithm works in three steps. First, a so-called 
energy envelope is calculated for each recorded waveform 
and a dynamic signal threshold is set. Then, for any 
signal with energy above this given threshold, a hit is 
recorded for the given antenna at that threshold-crossing 
time. In this way, hit patterns are generated for each 
event which, in the third step, are checked for consistency 
with incoming planar radio waves. 

The “energy” envelope is in principle calculated as the 
RMS of a sliding 5 ns time window of the voltage data 
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dt value/ a.u. 


(a) 



(b) 


Figure 16. (a) The hit time differences in histograms for each 
of the five geometrical groups for the event shown in Figure 
15(a) The pair groups of the histograms are schematically 


shown as inlay in the bottom plot. The quality parameter 
here (QP) = 1.6. (b) The hit time difference histograms for 
the noise event shown in Figure [l5(b)[ In this second case the 
quality parameter equals 0.5. 



0 5 10 15 20 25 

Time sequence QP // Hit count 


Figure 17. The quality parameter QP (solid line) compared 
to a simple count of all hits in a pattern as they appear in 
the example in Figure 15 (dashed line) for simulated neutrinos 
with energies between 10^® eV and 10^^ eV (blue) and thermal 
noise events (red). All values are scaled to cumulatively reach 
99% in both noise distributions at the same X value. The 
distributions are normalized to the total event count. 


separated in the presented algorithm. In total there are 
five groups of pairs with the same geometrical orienta¬ 
tion. For an incoming plane wave these histograms are 
expected to show a strong peak while they should be flat 
for thermal noise. The normalized sum of the maximum 
bin counts from each histogram is used as the time se¬ 
quence quality parameter (QP), to distinguish incoming 
wavefronts from thermal noise. 

The noi se rejection power of this filter is shown in 
Figure 17 for simulated neutrinos between 10^® eV and 
10^^ eV. In the range between 10^® eV and 10^® eV, 92% 
of neutrino signals are kept at 99.9% noise rejection. The 
actual cut has to be tightened to provide adequate ther¬ 
mal noise rejection for the full data sample. However, as 
will be shown, the signal efficiency remains high. 


B. Vertex reconstruction 


to enhance the signal-to-noise ratio. For each event and 
antenna channel a threshold is defined as /t^-l-4 • aE, with 
fiE being the average and aE the RMS of the full energy 
envelope of a waveform. Whenever an envelope crosses 
the threshold, a hit is recorded for the given channel with 
a coarse timing precision of 5 ns. The hits of all channels 
taken together form a hit pattern (Figure [T^. 

In the next step, pairs are formed from antennas at 
roughly the same depth (horizontal pairs) and anten¬ 
nas on the same string (vertical pairs). For pairs with 
the same geometrical orientation, the time difference, di¬ 
vided by the antenna distance, is filled into a common 
histogram as shown in Figure |16| This figure illustrates 
how radio signal patterns and thermal noise patterns are 


After separating impulsive radio signals from thermal 
noise, we use a directional vertex reconstruction to dis¬ 
tinguish neutrino-induced emission from anthropogenic 
noise. Here, we use the fact that man-made signals will 
reconstruct to the surface while only neutrino signals 
originate from within the ice itself. Furthermore, signals 
generated above the ice undergo refraction at the ice air 
boundary. Under the assumption of an index of refraction 
of 1.755 at the antennas, the critical angle at this inter¬ 
face can be calculated to be 55° and events produced in 
air are thus limited to a zenith angle between 90° and 55° 
as viewed in-ice from the ARA stations. Events which 
are generated directly at the ice boundary, for example 
by driving vehicles, will arrive with a minimal angle of 
roughly 40°. However, this is not a concern for winter 
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(b) 

Figure 18. (a) The cross-correlation graph for a calibration 
pulser signal using waveforms measured in two different an¬ 
tennas. (b) The two recorded waveforms after a shift by the 
time of maximum correlation dt = —78.4 ns. 


running, when on-ice activities are minimized. 

Background radio pulses from cosmic ray air showers 
can mainly be emitted from above the ice or as transi¬ 
tion radiation from particle bunches at the ice-air bound¬ 
ary [3H]. For the ARA stations they will appear in the 
same angular region as the anthropogenic background de¬ 
scribed above. In the background estimation they will be 
treated as surface events. Only penetrating high energy 
muon bundles might be able to generate dense enough 
cascades in the deep ice via catastrophic energy loss to 
produce detectable radio signals. Studies of UHE neu¬ 
trinos predict neutrino fluxes to be roughly 100 times 
smaller than the cosmic ray flux at 10^® eV m- How¬ 
ever, a study of the energy loss of muon bundles in South 
Pole ice below 1450 m indicates that at such a primary 
cosmic ray energy the probability for an energy loss of 
10^® eV over 5 m, the lower end of the ARA sensitivity, 
is less than 10“® [38]. This gives us confidence that the 
background expected from muon bundles is small com¬ 
pared to neutrino signals. The actual amount of back¬ 


ground, produced by such muon bundles in the ARA 
detector is currently under more detailed investigation. 

The reconstruction algorithm developed for this analy¬ 
sis is based on a system of linear equations, formed from 
the signal arrival times on the different antennas. In the 
ARA stations relative arrival time differences can be mea¬ 
sured using the cross-correlation g between two antenna 
waveforms fi and /2 (Figure [T^: 

9{t) = /i * /2 = mfi)r ■ Hf2)) , (10) 

where J- stands for the Fourier transform of a given func¬ 
tion. The maximum of the correlation graph should occur 
at the delay time between the two input signals. With 
this method, a timing precision of 100 ps on average can 
currently be achieved in ARA (see Appendix |A|) . 

In principle, the time differences from all possible pairs 
can be used in the reconstruction. However, to exclude 
antennas that did not register a signal waveform from 
the algorithm, the maximum correlation amplitude in the 
correlation graph is used as a selection criterion for good 
antenna pairs. For a pair to be selected, the squared 
correlation amplitude has to cross a dynamic threshold 
adapted to the overall signal amplitude in an event, but 
at least a fixed lower limit above the product of the inte¬ 
grated power of the two correlated waveforms. All pairs 
passing the threshold are initially used in the reconstruc¬ 
tion prior to a more refined channel pair selection per¬ 
formed in further steps, which is based on quality criteria 
applied to the outcome of the reconstruction. With the 
time differences found for good antenna pairs, a system 
of equations is set up, using the equality between the dis¬ 
tance to the signal source and the measured travel time: 

= {x„-Xif‘ (11) 

+ {vv - Vif 


where is the time of emission at the vertex, ti the 
time of reception by the antenna i and x,y and z the 
respective spatial coordinates. The speed of light c is 
assumed to be constant and equal to the average speed 
at the station depth c = 0.3/1.755 m/ns (from [25]). The 
changing index of refraction with depth is not taken into 
account in this reconstruction, which has an influence 
on the zenith reconstruction precision. This, as will be 
shown, does not have a significant effect on the efficiency 
of the analysis. 

When subtracting this relation for pairs of antennas, 
and after some reordering, one can obtain for each pair 
of antennas i and j: 


‘^Xij -j- y-i} • ‘^'Uij Zv ' ^ 


i;,ref 


2c^U 


= r/-r (12) 


Here, the index ij indicates the difference between the 
values of a certain parameter for antennas i and j, the 
parameter r denotes the distance to the center of the 
coordinate system and the index “ref” indicates a refer¬ 
ence antenna for which the signal arrival time is set to 
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Figure 19. The dependence of the azimuthal reconstruction 
on the residual (Equation 141. The Y-axis indicates the dif¬ 
ference between the reconstructed angle and the true angle 
for 165000 simulated neutrino events with energies between 
10^® eV and 10^^ eV. Events with a high residual triggered 
the detector but do not contain strong enough signal to be 
properly reconstructed and are likely thermal in origin. 


be to = 0. This relation is used to set up a system of 
equations, linear in the vertex coordinates and emission 
time, represented by the matrix equation 

Av = b, (13) 


with V containing the vertex coordinates and emission 
time and the matrix A and the vector b offsets and ar¬ 
rival time differences. It should be noted that this ap¬ 
proach is similar to Bancroft’s solution of GPS equations, 
described in [^. The solution of this equation can be 
obtained using matrix decomposition tools [iQl Hi] . It is 
thus not seed-dependent and very fast, which allows us 
to perform several thousand reconstructions per second. 
Their precision and stability depend strongly on the pre¬ 
cise knowledge of the relative antenna positions and pos¬ 
sible time offsets between their recorded signals, which 
can be caused by cables or other electronic components. 
A calibration of the station geometry and systematic time 
delays between antenna waveforms has been presented in 
Section [irC| For each reconstruction, a residual is calcu¬ 
lated as 


res 


b 

I 


A • V 


lA-i 


A, 


chp 


(14) 


with Achp the number of involved channel pairs. This 
indicates how well the reconstructed values fit the mea¬ 
sured arrival time differences. Based on the residual, the 
channel pair selection is refined to further exclude noise 
antenna pairs if they have not been identified in the first 
step of channel selection. 

The residual, alongside with some minor quality pa¬ 
rameters, indicates whether a reconstruction is consid¬ 
ered trustworthy. Figure 19 shows the dependence of 
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(a) 



(b) 

Figure 20. Simulated neutrino vertices reconstructed in az¬ 
imuth and zenith angle with all quality criteria applied, (a) 
The difference between reconstructed and true azimuth, (b) 
The reconstructed zenith angle plotted as a function of the 
true zenith angle of each event. 


the azimuthal reconstruction on the residual. For low 
residuals, a high quality reconstruction can be obtained, 
while reconstructions with a residual above 10“^'® ap¬ 
pear to point to locations that are broadly distributed 
in azimuth with respect to the true value. The residual 
also rejects more thermal noise events in favor of signal 
events although that is not its main purpose. 

The result of reconstructions from a set of simulated 
neutrinos after application of all reconstruction quality 
criteria is shown in Figure While for the azimuth 
reconstruction a precision of better than 2° (0.3° when 
excluding the tails) can be achieved, the zenith recon¬ 
struction is significantly degraded by surface reflections 

























15 


bO 

(U 

o 

u 

fV 


M 


-10 

-12 

-14 

-16 

-18 

-20 

-a 


ex 

pected /deg 

reco /deg 

ICl zenith 

- 22.7 

- 17.9 ±0.3 

IC22 zenith 

- 23.1 

- 18.3 ±0.3 

ICl azimuth 

260.0 

258.0 ± 1.3 

IC22 azimuth 

266.0 

264.0 ± 1.2 


l" 


T 


r 


250 


255 


260 


265 


Azim uth reco /deg 

(a) 


110 " 


10 " 


]\ 

270 


10" 


-10 

^_12 

"D 

O -14 

u 

(V 

-F -16 
-18 

-20 

-22 


d) 

M 


ex 

pected /deg 

reco /deg 

- ICl zenith 

- 19.5 

- 15.1 ±0.3 

IC22 zenith 

-20.6 

- 16.0 ±0.4 

ICl azimuth 

230.0 

231.0 ±0.2 

IC22 azimuth 

234.0 

235.0 ±0.3 




, 10 " 


10 " 


226 228 230 232 234 236 238 240 

Azim Uth reco /deg 


i 


10 " 


(b) 


Figure 21. Reconstruction of pulsers deployed in IceCube 
holes in the deep ice (a) with A2, (b) with A3. The plots 
include data for the expected positions and the reconstructed 
positions with their standard deviation. The influence of sys¬ 
tematic differences of a few degrees between the true and re¬ 
constructed angles on neutrino identification is negligible. 


and ray-tracing effects. About 30% of the events show 
behavior which causes the reconstruction to miss by sev¬ 
eral degrees in zenith angle. This causes an efficiency 
loss of roughly 6% due to application of our angular cuts 
on radio waves coming from the surface, which is not 
dramatic. 

Figure [TT] shows the reconstruction of two pulsers de¬ 
ployed at a depth of 1450 m in the ice on IceCube strings 
1 and 22 deployed in the final season of IceCube con¬ 
struction. Their distance to the ARA stations is roughly 
4000 m. These pulsers are therefore our most neutrino¬ 
like calibration tool. The plots show that both pulsers 
can be reconstructed with good precision by both sta¬ 
tions. In addition to the roof pulser, this is another 


external source which confirms that the reconstruction 
algorithms should work properly for neutrinos. 


C. Cuts and background estimation 

Based on the presented algorithms, three cuts are used 
to distinguish neutrino signals from thermal and anthro¬ 
pogenic background. Thermal noise must be reduced by 
a factor 10“^° to reach the goal of ten times less events 
than the expected number of neutrinos. This can be 
achieved mainly by requiring a time sequence quality 
parameter of at least 0.6 to select an event. Further- 



fa) 



Figure 22. (a) The distribution of recorded data in the two 
main cut parameters. All data correlated to known radio 
source locations (calibration pulser, surface) by reconstruc¬ 
tion are removed. Hence, the remaining events can be ther¬ 
mal noise, mis-reconstructed radio events or neutrinos, (b) 
The distribution of simulated signal (see Section III I in the 
two main noise parameters. 
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more, to select high quality reconstructions and reject 
the remaining thermal noise, the residual is required to 
be less than 10“"^. In this way only well-reconstructed 
impulsive radio signals are kept, which can further be re¬ 
duced by angular cuts. As the cut values for the time 
sequence parameter and the reconstruction residual have 
not been optimized in a strictly systematic way, we note 
that there might be room for improvement in a subse¬ 
quent analysis employing these algorithms. Figure 
illustrates the steep decline of the noise distribution to¬ 
wards the cut boundaries in the two main cut parame¬ 
ters. As a representative thermal noise sample, all data 
from the 10% burn-sample which cannot be correlated 
to a known signal source (pulsers, surface activities) are 
used. Simulated signal events from the dataset described 
in Section |III| are distributed broadly compared to that 
noise sample. 


Angular cuts are placed around the known locations 
of calibration pulsers inside the ice and are specific to a 
given station. In addition, a surface cut is applied, reject¬ 
ing all events reconstructed to a zenith angle oi 9 > 35° 
for A2 and 9 > 40° for A3. This cut can be a bit looser 
for A3 since the reconstruction errors are smaller for this 
detector (see Figure 21). The reason for this difference 
in precision is due to the different number of channels 
in both stations which are available for reconstruction. 
Whereas all 16 channels are used in station A3, in A2 
only 13 channels are included in vertex reconstruction. 
One channel, D4BH, is broken, while two other chan¬ 
nels show a puzzling timing offset which could not be 
removed in the geometrical calibration. The cut values 
are chosen for each angular requirement separately, such 
that each allows less than 0.01 background events to en¬ 
ter the signal sample in the full data set. The num¬ 
ber of background events expected to pass a given cut 
is estimated from the 10% data subset by fitting an ad¬ 
equate Gaussian or exponential function to the tail of 
an event distribution close to each cut. The best ht pa¬ 
rameters and the position of the cut are used to obtain 
the number of background events expected to leak from 
the angular region being excluded. The uncertainty on 
the number of background events is derived from the fit 
errors. Note that calibration pulser events are normally 
tagged by the DAQ as calibration events and excluded 
from the analyzed data sample. However, due to pos¬ 
sible mis-tagging, pulser events may leak into the final 
sample. Therefore, all pulser events are taken into ac¬ 
count in the background estimation, even if they are not 
part of the 10% burn-sample. This is a very conservative 
estimate but strengthens the analysis against mis-tagged 
pulser events. 


After these cuts, the background expectation for the 
full data recorded in the year 2013 is 0.009 ± 0.010 for 
A2 and 0.011 ± 0.015 for A3. The number of neutrinos 
expected to be observed by the combined two-station de¬ 
tector from the flux prediction in m for a crossover en¬ 
ergy from Galactic to extragalactic cosmic ray sources of 
Emin = 10^®-® eV amounts to 0.10 ± 0.002{stat) events. 
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Figure 23. Reconstructed events that passed the thermal 
noise and reconstruction quality cuts for (a) A2 and (b) A3. 
The black boxes indicate the angular cut regions around the 
calibration pulser positions and the black line indicates the 
surface cut. Events inside the squares and above the surface 
line are rejected. 


V. RESULTS AND CROSS CHECKS 
A. Results 

The results of the above described analysis are sum¬ 
marized in the two sky maps in Figure No events are 
found outside the angular cut regions which implies that 
no neutrino candidates have been observed. This agrees 
with the expectation of 0.1 signal and roughly 0.02 back¬ 
ground events in the two stations. The difference between 
summer and winter source locations at South Pole is pre¬ 
sented in Figure |2^ showing the impact of human activ¬ 
ities during summer which is limited to surface events. 
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Figure 24. Reconstructed events that passed the thermal 
noise and reconstruction quality cuts for A3 (a) in the austral 
winter and (b) in the austral summer when station activity 
is maximal. Note that due to improper tagging of calibration 
pulser events in the first months of winter, many such events 
entered the final data sample. This is however accounted for 
in the final angular cuts of the analysis. 


With this result from two ARA stations (designated 
as “ARA2”) we can calculate a differential limit on the 
neutrino flux in the sensitive energy region as shown 
in Figure |25[ For the neutrino energy range between 
10^® eV and 10^® eV, the energy where most neutrinos 
are expected to be observed, the limit is E^Fup{E) = 
1.5 X 10"® GeV/cmVs/sr. 

This limit is calculated as 


E^E^piE) 


E^S{E) 
ES{E) ■ 


K{E) 

dE 

K{E) 

In(lO)’ 


(15) 


where the factor K(E) is derived with the construction 
described in m as the 90% Poisson confidence limit for 
no observed events under the expectation of zero back- 



Figure 25. Neutrino limits and sensitivities from various ex¬ 
periments including the 7.5 months data analysis of the two 
ARA stations described herein. Systematic errors, as derived 
in Section lYil have been accounted for in the ARA2 limit. 
The ARA37 (3 yr) sensitivity is projected from the ARA2 
(7.5 m) trigger level without accounting for systematic errors 
{K{E) = 2.44). Data for other experiments are taken from 
[71 [^l42ll45] . The neutrino fluxes are derived in [37] and [46| . 
See Appendix [B] for EF{E) scaling on the y-axis. 


ground. Systematic errors on the signal efficiency, as 
described in Section [VTl are accounted for in this factor 
following the method presented in [48j with the improve¬ 
ment proposed in |49j . The error caused by uncertainties 
on the cross section is not taken into account. S{E) de¬ 
notes the sensitivity of the detector which is calculated 
from the effective areas Aqs and live times T of each 
detector as 

= WW4— 

47r • (Aeff,2 • T 2 + ^eff.3 ' Ts) 

Furthermore, the limit is presented as a half decade inter¬ 
polation for a logarithmic energy scale with a resolution 
of dLogio(if) = 1 logarithmic bins. Therefore we obtain: 

dE = E- ln(10)dLogio(Fl) = E ■ In(lO). (17) 

The resulting limit for two ARA stations 10 months after 
deployment is not yet competitive with the current best 
limits from the IceCube detector. In spite of this they 
show, when projected to the full size of ARA37, that the 
completed detector is expected to be sensitive to main¬ 
stream models for neutrinos from the GZK process. 


B. Cross checks 

Although the limit obtained with the currently de¬ 
ployed ARA stations is not competitive yet, this first data 
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Figure 26. A South Pole map with the two ARA stations and 
the XY reconstruction of an A2/A3 coincident event series via 
the parallax method (green dots). The events in this series 
are distributed over a time of 50 s starting from the rightmost 
point in the map. 


analysis proves the capabilities of the full ARA detec¬ 
tor. Cross checks have been performed with events that 
are observed coincidentally in both stations to demon¬ 
strate that the employed algorithms select impulsive ra¬ 
dio signals and that the directional reconstruction works. 
Events that have passed both the thermal noise and re¬ 
construction quality cuts are considered coincident if they 
trigger both stations within a time window of 12 ps. This 
time corresponds to the maximal in-ice travel time of sig¬ 
nals between the two stations given their separation of 
2 km. Sequences of such events, appearing in short time 
frames of roughly 60 s have been found which show evi¬ 
dence of originating from airplane communication trans¬ 
mitters. 

With the azimuthal reconstruction from each station, 
the XY position of the source can be determined via 
trigonometric calculations. The crossing point of the two 
beams pointing to the azimuthal reconstruction of both 
stations is used as its XY position. This method is in 
the following referred to as the “parallax reconstruction”. 
One particularly interesting event sequence is shown in 
Figure The positions of the events within this hit 
series form a smooth track, indicating an emitting object 
that moves at a speed of several hundred km/h at an 
increasing height of order 500 m above the ice surface. 
Waveforms recorded for events at different points of the 
sequence are presented in Appendix This track is very 
useful as a cross check since it passes on top of station 
A2 which should be evident in the zenith reconstruction. 
In Figure!^ the expected zenith angle from the XY po¬ 
sition is compared to the zenith reconstruction of station 
A2, showing good agreement within the error bars. 


Figure 27. Comparison between the expected zenith from the 
A2/A3 XY parallax reconstruction (green dashed line) and 
the reconstructed zenith angle by A2 only (green solid line). 
Errors from the XY reconstruction, calculated from the errors 
on the reconstruction of a surface pulser with known position, 
are shown as a green band. 


The track confirms that the detector is capable of ob¬ 
serving radio sources. Furthermore, the agreement be¬ 
tween different positioning methods and the smoothness 
of the reconstructed track are evidence that the used 
analysis tools work properly to identify such sources and 
reconstruct their position. 


VI. SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties of the presented analysis 
result from errors on the theoretical models of neutrino 
interactions and radio wave propagation, as well as on 
the calibration of the detector. The error estimation is 
performed in a similar way as has been described in [26] . 

The neutrino interaction cross section at the energies 
of interest above 10^® eV is calculated based on measure¬ 
ments at much lower energies and is thus subject to large 
uncertainties. To check the influence of this uncertainty 
on the effective area of the ARA detector, simulations are 
run using the upper and lower limits of the cross section 
estimates from [33]. The effect on the final analysis re¬ 
sult can be seen in Figure]^ Especially for the highest 
energies, this is the dominant uncertainty in the analysis. 

A further uncertainty results from the error on the ra¬ 
dio attenuation length measurement in the South Pole 
ice sheet. This measurement has been performed using 
calibrated pulsers, deployed at different depths with the 
last IceCube strings [21]. With the obtained data, the 
local attenuation length at a given depth can be inferred 
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Figure 28. The relative difference in effective area at analysis 
level, caused by various systematic error sources. 


with knowledge of the temperature and density profile of 
the ice. The difference between this result and an earlier 
measurement, using the bed rock under the ice as a re¬ 
flector for radio waves emitted and received at the surface 
[50], provides a measure of the uncertainty in the atten¬ 
uation length. The error on the effective area is again 
obtained by comparing simulations with different sets of 
parameters. As visible in Figure 28 it contributes only 
slightly to the final error. 


One should note that the uncertainty on the changing 
index of refraction inside the ice is not a major concern 
for ARA due to the deep deployment of the stations at 
180 m below the ice surface. Below this depth the index 
of refraction does not change appreciably and approaches 
the value of 1.78 for the deep ice. 


One more important uncertainty relates to the final 
signal to noise ratio recorded by the signal chain. This 
depends on the assumed ambient noise, the antenna 
directivity, the transmission coefficient and the assumed 
noise figure as explained is Section |II D 3[ To estimate 
the resulting systematic error on the effective area, 
the overall amplitude of the incoming radio signal is 
reduced according to the results from Section |II D 3| 
Under these conditions, sets of neutrino events are 
simulated and the analysis is re-run. The presented 
analysis is nearly exclusively based on coarse envelopes 
of time domain waveforms or time differences, derived 
from cross-correlation. Therefore, precise knowledge 
of the frequency response is of secondary importance. 
Figure shows that the uncertainty on the signal chain 
calibration has the greatest impact at low energies, 
when most of the incoming signals are weak and the 
signal to noise ratio is low. At higher energies, this error 
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Figure 29. The residual of reconstructed events as function of 
the SNR for examples of (a) neutrino simulations in A3 and 
(b) calibration pulser signals recorded by A3. 


loses importance compared to the error on the cross 
section. Due to the limited knowledge about the upper 
systematic errors in that measurement, as explained 
in Section [IIP 3[ such a limit is currently not quoted. 
Only the lower error is used in the determination of the 
neutrino limit presented in Figure ]^ 

The last estimated uncertainty is the difference in the 
analysis efficiency obtained from simulated neutrinos and 
recorded data. The shape of noise and signal waveforms 
will not perfectly match between simulations and real 
data. Therefore, the used analysis algorithms are com¬ 
pared between real data from calibration sources and 
simulated neutrino signals. For this comparison, the sig¬ 
nal to noise ratio (SNR) is chosen as a simple parameter, 
independent of the shape of the waveform. In Figure 
calibration data from station A3 are compared to 
simulated neutrino events with energies between 10^® eV 
and 10^^ eV for one of the two main cut parameters, the 
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Figure 30. Top panel: The analysis efficiency of the ARA 
stations for simulated events (green) and for calibration pulser 
events on A3 (blue). Bottom panel: The relative difference 
in efficiency between simulation and pulser data (red) and 
the normalized event count Ni per SNR bin (orange). The 
shown simulated events are produced with a primary energy 
of 10^® eV. 


reconstruction residual. Generally, the distribution for 
simulations aligns well with the recorded data. The dif¬ 
ference around the cut value seems especially small. 

For the systematic uncertainties it is interesting to 
quantify the analysis efficiency, i. e. the fraction of events 
that pass the applied cuts, for real data and simulations 
and to calculate their relative difference. To do so, one 
can compare its dependence on SNR after separating the 
simulated neutrinos into single decade energy bins. Fig¬ 
ure!^ shows this comparison for a neutrino energy of 
10^® eV. The difference |Acff,i| is plotted relative to the 
efficiency in simulation. For each energy an average dif¬ 
ference in efficiency, weighted by the number of events in 
each SNR bin, is calculated as 


a 


sys — 


(l^eff,»| ■ N^) 

ENi 


(18) 


with Ni the normalized number of events per SNR bin. 
This difference is propagated as a systematic error on the 
analysis efficiency into the limit calculation. As visible 
in Figure this is a non-negligible contribution to the 
systematic error. However, with a better knowledge of 
the detector, simulations will become more precise and 
this uncertainty can be reduced. 

The systematic uncertainties calculated for the number 
of recorded neutrinos are summarized in Table m 


Table II. The systematic error from various sources on our 
signal expectation of 0.10 neutrino events for a flux prediction 
from [37]. 


Source 

positive error negative error 

Cross section 

0.035(34%) 

0.020(19%) 

Attenuation 

- 

0.005(5%) 

Signal Chain 

- 

0.011(10%) 

Analysis efficiency 

0.018(17%) 

0.018(17%) 

TOTAL 

0.040(38%) 

0.030(29%) 


VII. SUMMARY AND OUTLOOK 


We have demonstrated the power of ARA as an UHE 
neutrino detector through an analysis of data from the 
two deep stations currently in operation. Through cal¬ 
ibrations, good timing precision and geometrical under¬ 
standing of the detector have been achieved, which are 
key factors in the detection of radio vertices. Further¬ 
more, initial analysis algorithms have been presented. 
These show a good efficiency for retaining signal (60% 
efficiency at 10^® eV) with background rejection by 10 
orders of magnitude, leaving 0.02 events, at the current 
trigger settings. Thermal noise can be rejected by sim¬ 
ple algorithms to a high level and radio vertices can be 
reconstructed with an angular precision of a few degrees. 
This relatively simple reconstruction algorithm is found 
to have a RMS of < 2° in azimuth reconstruction and 
appears to be stable for most cases of zenith reconstruc¬ 
tion (see Figure 201 without accounting for ray tracing 
effects. Improvements and alternatives to the algorithm 
are currently being developed. 

In addition, cross checks confirm that the used analy¬ 
sis algorithms select radio signals from background and 
return sensible directional reconstructions. 

Since the probability for a neutrino detection within 
the given time of operation was very low and of order 
0.10 events, the analysis algorithms have not been op¬ 
timized for reconstruction of neutrino four-vectors. To 
determine the incoming direction and energy of a neu¬ 
trino, one needs to know the polarization of the event 
and the distance to the vertex. The polarization can be 
determined by comparing the signal strength in Hpol and 
Vpol. The distance reconstruction is important for the 
determination of the neutrino energy. Since the distance 
reconstruction is not crucial for neutrino identification, 
it is not essential to the ARA analysis strategy. For now, 
only a lower limit on the distance can be provided for 
each event, which results in a lower limit for the neutrino 
energy. The achievable energy resolution of the ARA 
detector is thus still unknown and more work in this di¬ 
rection is needed to prepare for the case of a neutrino 
detection. 

The presented limit, resulting from one year of data 
taking, is not yet significant but raises the expectations 
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fore expected to show a significant improvement in 
precision and reconstruction efficiency. 

• The detector calibration is under continuous im¬ 
provement which will help especially to reduce sys¬ 
tematic uncertainties on the detector geometry, 
timing and the signal chain, which consequently 
improve analysis algorithms like the angular recon¬ 
struction. 


Figure 31 shows the numbers of neutrinos which are 
expected to be seen with the ARA37 detector, given the 
trigger and analysis efficiency for 2013, within 3 years of 
operation at different energies and from different fluxes. 
Additionally, an expectation for a power law flux, nor¬ 
malized to the Waxman-Bahcall bound from m is plot¬ 
ted to illustrate the response of the detector over a wide 
range of energy. This figure shows that the planned full 
ARA detector is a promising candidate for the detection 
of ultra-high energy neutrinos. 


Figure 31. Numbers of neutrinos versus energy in half-decade 
bins, as expected to be seen with the ARA37 detector within 
3 years. These numbers have been extrapolated from the 
presented trigger and analysis efficiency for the first two sta¬ 
tions. The numbers are calculated for 3 different ffux predic¬ 
tions which are partly shown in Figure and estimated in 
|37[ I46[ I51| . The total expected number of neutrinos for the 
predictions are 9.4 (Ahlers 2010) and 5.5 (Kotera 2010). 

for the discovery of UHE neutrinos with the full ARA 
detector. 

The results presented in this paper reflect the sta¬ 
tus of detector operations in 2013 and of currently 
available analysis tools. Since then, several improve¬ 
ments have been developed, benefiting coming analysis 
results, of which the most important are the following: 

• The live time per year has been increased by 
roughly 25%, thanks to newly developed monitor¬ 
ing tools which allow for quick debugging, particu¬ 
larly in the case of downtime due to issues in the 
detector electronics. 

• The trigger and readout windows have been opti¬ 
mized and widened to enhance the detector sensi¬ 
tivity and to render analysis tools, complicated by 
cutoff signal waveforms, more efficient. 

• A PCI Express bus has been integrated into the 
DAQ to replace the previously used USB connec¬ 
tion. This allows the recording of data at an event 
rate several times higher than before, and allows 
lowering the trigger threshold, enhancing the neu¬ 
trino sensitivity. 

• New reconstruction algorithms are underway which 
take ray-tracing effects into account and are there- 
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Appendix A: Calibration of the IRS2 digitizer chip 

The 1RS2 chip is a custom ASIG for radio frequency ap¬ 
plications [S2] . It is designed to digitize at a speed of sev¬ 
eral GS/s at low power consumption of less than 20 mW 
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per channel. These features merit usage of CMOS tech¬ 
nology for the implementation of the sampling and digiti¬ 
zation steps, as well as utilization of a large analog buffer, 
which allows for a slow digitization technique without in¬ 
ducing dead time. In the IRS2 chip the data are sam¬ 
pled via a Switched Capacitor Array (SCA) which utilizes 
finely tuned delay elements to set up a sampling sequence 
of the input data. These delay elements can differ from 
their nominal delay width due to variations during the 
chip fabrication and have to be calibrated individually. 
The SCA consists of 128 sampling capacitors per channel, 
equally divided into even and odd samples on two delay 
lines and each with a delay element requiring individual 
calibration in timing. In addition to that the ADC to 
voltage conversion gain needs to be determined for each 
of the 32768 buffer elements on each channel to obtain a 
proper voltage calibration. 


For both calibrations, timing and ADC gain, sine 
waves are recorded and fit to ideal waveforms (Figure 


32). The data are recorded in the laboratory with the 


instrument at the final in-ice temperature of —50 °C. 


For the timing, first the fit frequency is compared to 
the frequency of the input waveform to calibrate the av¬ 
erage sampling speed. Then, each individual sample tim¬ 
ing is compared to the fit waveform to obtain a correction 
factor for the given delay element. This correction is ap¬ 
plied if samples have an absolute ADC count below 30, 
thus if the derivative of the sine wave is maximal and 
the influence of voltage errors is small. Corrections are 
directly applied and the process is repeated for several 
iterations until the correction factors converge. Figure 
1^ shows the final timing corrections needed for samples 
of a selected channel after several iterations of calibra¬ 
tion. One can observe a statistical fluctuation symmetric 
around zero, which indicates that all systematic errors 
have been removed by the calibration. The underlying 
calibration values are then used to set the correct timing 
for each delay element. The visible spread of some dis- 



Figure 32. A typical calibration waveform separated by even 
(blue) and odd (green) samples with a fit waveform (red). 
The horizontal lines indicate the range for individual sample 
correction in timing. 



(a) Even samples 
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Figure 33. The timing corrections for all delay elements (X- 
axis) of one channel in station A3, separated by an even (a) 
and odd (b) delay line after several iterations of calibration. 


tributions is connected to a non-linearity in the voltage 
of the respective sample. 

For the voltage calibration, timing-corrected wave¬ 
forms are used as input. Due to the density of the chip 
structure, a slight non-linearity and asymmetry around 
0 is induced in voltage which depends strongly on the 
channel number. Furthermore, we need calibration data 
for about 650000 storage elements per station, which re¬ 
quires a huge data sample. These conditions render clas¬ 
sical voltage calibration methods difficult. For the cal¬ 
ibration of the ADC to voltage conversion gain of the 
IRS2, input sine waves of known amplitude and frequency 
are fit and ADC samples are compared to the fit whenever 
its derivative is smaller than 45% of the maximal value. 
In contrast to the timing calibration, a small derivative 
is required, this time to minimize the influence of timing 
errors. Following this procedure a statistically significant 
sample for each storage element can be collected, using 
input waveforms over a wide range in amplitude (Figure 

Further calibrations which have been performed on the 
IRS 2 chip are a check of 

• the frequency response, which could not be deter- 
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Figure 34. Calibration data for the calculation of the ADC- 
to-voltage conversion for a single storage sample: Collected 
data (red dots), averaged data (black points) with errors and 
a broken 3rd order polynomial fit to the average data (blue 
line). 



time delay / ns 


Figure 35. Arrival time difference for a calibration pulser 
signal on two channels of station A2 after timing calibration 
(red) and the full calibration (blue). 


mined conclusively with the used data set, 

• the temperature dependence of timing and voltage, 
which appears to be negligible in the temperature 
range of the ARA experiment. 


The main purpose of the calibration is to obtain good 
correlation timing between incoming waveforms. There¬ 
fore, the calibration is cross checked with calibration 
pulser waveforms, recorded on different channels. On 
average, a precision of 100 ps can be achieved (Figure 


35), which is entirely adequate for good vertex angular 


reconstructions. Determination of the radius of curva¬ 
ture of the incoming wavefront is considerably more dif¬ 
ficult with stations of limited size; therefore, for sources 
more than tens of meters from the station, the range to 
emission is effectively an unknown at the current level of 
analysis. 

The influence of the timing jitter and possible non¬ 


linearities in the ADC to voltage conversion gain are 
tested on a simulated waveform, sampled at 3.2 GS/s, 
with a frequency spectrum similar to what is expected 
for an ARA signal. Timing calibration errors are modeled 
as Gaussian distributed random jitter for different stan¬ 
dard deviations and added to the sample timing. Non- 
linearities, left in the calibration of the ADC to voltage 
conversion gain, are modeled with a simple third order 
polynomial as 

gdist = g +k-{g‘^ + (Al) 

with g the linear relation between the ADC count and 
the input voltage and k the level of a non-linear addi¬ 
tion. As for the timing jitter, k is a, random number for 




(b) 

Figure 36. Comparison of Fourier spectra for a simulated 
waveform (undisturbed in blue) at various error levels in the 
digitizer calibration, (a) Different magnitudes of timing jitter 
applied to each sample. The RMS represents the width of the 
Gaussian distribution used to generate the random jitter for 
each sample, (b) Different levels of a non-linear addition in 
the ADC to voltage conversion gain as presented in Equation 
|A1| The RMS represents the width of the distribution used 
to generate k. 
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each sample following a Gaussian distribution of a given 
width. The resulting Fourier spectra of the original and 
smeared waveform for different error levels are presented 
in Figure As is visible, the smearing adds a broad¬ 
band component to the spectrum which has a small in¬ 
fluence on regions of strong signal but adds significantly 
at low signal levels. The shown levels of smearing have 
been chosen to be within a reasonable range as observed 
during the calibration process. Additional investigations 
are needed to determine the actual amount of jitter and 
non-linear gain components for each sample to quantify 
the influence on the Fourier spectrum more precisely. 

Appendix B: Neutrino limits with alternative EF{E) 
scaling 

Figure shows the same data as Figure with the 
Y-axis scaling changed from E‘^F[E) to EE{E). This 
plot has been added for convenience of the reader. 



Figure 37. The neutrino limits and fluxes of Figure with 
an alternative y-axis scaling. 


Appendix C: A2 events reconstructed to different zenith angles 


In the event sequence for which sources are observed coincidently in both stations, as explained in Section [VB[ the 
zenith angle towards A2 changes drastically over the development of the track. This allows for observation of a shift 
in polarization in these events. Figures [38l |39| and |40| show a sequence of three events drawn from three different parts 
of the track. As expected, the event emitted vertically above station A2 shows a strong Hpol component, which is 
considerably weaker in those two events observed at shallower incident zenith angles. 
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Figure 38. Recorded waveforms for an A2 event, which was part of the event sequence presented in Figure [M] The event shown 
corresponds specifically to the source at the point indicated by the red X (bottom right plot), as it moves above the array. The 
reconstructed zenith angle in A2 is 51.6°. Note that channel D4BH on A2 is not operational. 
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Figure 39. Recorded waveforms for an A2 event, which was part of the event sequence presented in FigureThe event shown 
corresponds specifically to the source at the point indicated by the red X (bottom right plot), as it moves above the array. The 
reconstructed zenith angle in A2 is 87.0°. 
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Figure 40. Recorded waveforms for an A2 event, which was part of the event sequence presented in Figure [M] The event shown 
corresponds specifically to the source at the point indicated by the red X (bottom right plot), as it moves above the array. The 
reconstructed zenith angle in A2 is 57.5°. 
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